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The electronic structure of nanoscopic oxide-coated aluminum islands is investigated using a tight- 
binding model that incorporates the geometry, chemistry and disorder of the particle. The oxide 
coat is found to significantly increase the volume accessible to electrons at the Fermi level. The level 
statistics agree with random matrix theory predictions. States near the Fermi level show pronounced 
clustering regardless of disorder. It is suggested that the observed clusters of tunneling resonances 
may have a more complex origin than if they were solely due to many-body non-equilibrium effects. 

PACS numbers: 73.22.-f, 73.22.Dj 
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It is expected that the statistics of the discrete energy 
spectra of small disordered metal grains should be de- 
scribed by random matrix theory (RMT) jjj . Experimen- 
tal confirmation of RMT predictions for such systems, 
however, presents a challenge. Recently, some of the 
technological difficulties have been overcome and elec- 
tron tunneling spectroscopy experiments have been car- 
ried out on Al, Au and Co nanoscopic particles 
These transport experiments have opened the possibility 
of testing directly the applicability of RMT to the energy 
spectra of these nanoscopic metal islands (nanoislands). 

The observed electron tunneling spectra of Al nanois- 
lands coated with aluminum oxide |]2| showed a surpris- 
ingly high density and pronounced clustering of conduc- 
tance resonances. Agam and co-workers || proposed that 
this unexpected behavior is a manifestation of the ef- 
fects of electron-electron interactions on non-equilibrium 
states of the nanoislands. By focusing on the many-body 
aspects of the problem, they were able to account for 
the anomalous features of the data qualitatively in an 
appealing way. However, they modeled the underlying 
one-electron spectra of the nanoislands phenomenologi- 
cally. For a complete understanding of the experiments a 
microscopic treatment of the single-electron energy spec- 
tra is desirable. Microscopic models of aluminum nano- 
structures have been proposed IJ, but none of those 
works addressed the electronic structure or transport 
properties of Al nanoislands similar to those in the above 
experiments. Recently, Campbell et al. performed state- 
of-the-art molecular dynamics simulations of the passiva- 
tion of Al nanosphcres with aluminum oxide but did not 
calculate the electronic spectra of the particles |t| . 

In this communication, we present microscopic tight- 
binding calculations of the single-electron electronic spec- 
tra of Al nanoislands coated with Al-oxide. We include 
the geometry, atomic lattice, chemistry and disorder of 
the nanoparticle explicitly in our theory. Our leitmotiv is 
twofold: i) comparison of RMT predictions for the energy 
spectra of generic disordered metal grains with the elec- 
tronic structure obtained from realistic numerical simu- 
lations; and ii) improvement of our understanding of the 
experimentally observed tunneling resonance spectra. 



Motivated by the experiment, we consider a hemi- 
spheric nanoparticle of volume V ~ 40nm 3 , and investi- 
gate the effects of the disorder at the metal/oxide inter- 
face on the electronic structure and transport properties. 
(Different geometries yield qualitatively similar results.) 
We start with an analysis of the mean level spacing and 
level fluctuations in the neighborhood of the Fermi en- 
ergy of the particles. We find that describing the discrete 
energy levels in this region according to RMT should be 
a satisfactory starting point when dealing with Al nanois- 
lands coated with Al-oxide. However, the magnitude of 
the calculated mean level spacing is consistent with a sig- 
nificant enhancement of the effective electronic volume 
of metallic nanoparticle due to the presence of the ox- 
ide coat. Next, a closer look at a few levels around the 
Fermi level reveals a rich electronic structure: There are 
clusters of levels and the clustering is strongly affected 
by the specifics of disorder. Furthermore, the clustering 
of the single-electron energy levels persists even in the 
presence of severe interface disorder. Finally, we model 
the differential conductance (dl/dV) spectrum and find 
clusters of tunneling resonances. The clustering is due to 
the complex structure of energy levels around the Fermi 
energy of the particle. The clusters of resonances are 
sensitive to disorder and to the capacitances of the tun- 
nel junctions. The predicted number of tunneling reso- 
nances is not sufficient to explain the experimental data 
H within a purely single-electron picture of transport. 
Nonetheless, our results suggest that the observed clus- 
ters of tunneling resonances may have a more com- 
plex origin than if they were solely due to many-body 
non-equilibrium effects during transport. 

We assume that the Al atoms in the interior of the 
nanoparticle form an fee lattice with a lattice parameter 
a = 0.405nm. The atoms at the surface of the nanopar- 
ticle have a different coordination number from the bulk 
fee lattice. This establishes a criterion for defining the 
surface of the particle. In addition, the Al nanoislands in 
the experiment were passivated with a thin Al-oxide layer 
which acts as a tunnel barrier between the particle and 
leads. This layer also imposes some structural disorder at 
the metal-oxide interface. We model this metal-oxide in- 
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terface in the following way: First, we define which atoms 
of the particle belong to the surface and identify this re- 
gion as the interface; second, we randomly choose some 
of those sites (approximately 50%) to represent O atoms 
while the others correspond to Al. In the first part of 
this procedure we adopt two different criteria: Criterion 
1 (Ci): We assume the surface atoms to be those whose 
first neighbor coordination number is not equal to that of 
an fee lattice. Criterion 2 (C 2 ): Surface atoms are those 
whose first and/or second neighbor coordination num- 
bers are not equal to those of an fee lattice. The random 
choice of O and Al atoms in the metal-oxide interface is 
the only source of disorder in the results presented here. 

We use a multiband (s, p, and d valence orbitals) tight- 
binding Hamiltonian without spin-orbit coupling: 

H = ZiAcfra + T i a M { C h c j<* + 4^) w 

i,j label the atoms, c\ a (ci a ) creates(destroys) an elec- 
tron on site i, and the index a((3) indicates the s, p or d 
orbital, ef and Tjajp are the Slater-Koster (SK) || on- 
site and hopping (up to second neighbors) parameters. 
We adopt an SK model in which the atomic orbitals on 
different sites are orthonormal. For the Al sites in the 
interior (bulk) of the particle we used the SK parameters 
given in Ref. H. The Al and O sites in the metal-oxide 
interface have different SK parameters as is discussed be- 
low. The bonding of the atoms that form the oxide layer 
is due to charge transfer from the Al to the 0. Camp- 
bell et al. have shown that the metal-oxide interface of 
spherical Al nanoparticles consists mainly of intercalated 
O -1 / 2 and Al +1 / 2 j7|. These valence charges determine 
the SK on-site parameters for the interface atoms which 
we calculate following the Mulliken-Wolfsberg-Hclmholz 
molecular-orbital approach . We determine the near- 
est neighbor O-Al hopping parameters by assuming an 
average separation of 0.18nm between O and Al atoms 
in the oxide, applying Harrison's model to obtain the 
two-center transfer integrals and then transforming 
them to find the SK hopping parameters |9). The same 
procedure is applied to find the nearest neighbor 0-0 
hopping, using 0.30nm [l2]| as the average separation be- 
tween atoms. The SK hopping parameters between 
Al atoms are assumed to be the same regardless of the 
valence charge. 

We now present our results for a nanoisland of vol- 
ume V = 40.45nm 3 for different realizations of disor- 
der and surface criteria Ci and C2. For this volume 
the particle-in-a-box (SP) mean level spacing 8g l P = 
(4Ef /3A0V- 1 = 2.09meV where Ef is the Fermi en- 
ergy and J\f the electron density of bulk Al. The total 
number of atoms is N — 2587 and the surface-to-bulk 
ratio is (N s /N B )i = 0.538 and {N S /N B )2 = 0.970 for 
criteria Ci and C2, respectively. Note that for Ci the 
nanoisland is primarily bulk while for C2 it is almost 
evenly balanced between bulk and surface. In the lat- 
ter case the disorder is particularly severe. In Table | we 



present, for both Ci and C2, the number of O -1 / 2 (No) 
and Al +1 / 2 (Nai) atoms, the total number of electrons 
(n) in the nanoisland, the calculated Fermi energy (Ef), 
and the calculated mean energy level separation (8). Ep 
is the energy of the highest occupied level assuming that 
spin up and down levels are degenerate. We find that the 
value of Ep does not show major changes when switching 
from Ci to C2. It seems that the increase in n, due to the 
increased metal-oxide interface, in going from Ci to C 2 is 
compensated by the associated changes in the electronic 
structure of the particles. This compensation keeps Ep 
almost independent of the surface criteria. Surprisingly, 
we see that the different choices of surface criteria have 
a minimal effect on the values of S. Additionally, these 
values of 5 agree well with 8g l P . These two findings are 
non-trivial and important: They imply that the mean 
level separation around Ep is determined not by the vol- 
ume of the aluminum core of the nanoparticle but by a 
larger volume that includes part of the oxide coat. For 
small particles with large surface to volume ratios this 
should have a significant effect (exceeding a factor of 2) 
on the mean level spacing which should be reflected in 
the experimental spectra. However, a qualitatively sim- 
ilar enhancement of the value of the nanoisland volume 
that is deduced from the experimental capacitance mea- 
surements of the tunnel junctions may also be possi- 
ble. In this work the total volume V that determines 8g l P 
is kept fixed at an experimentally estimated value for the 
different disorder realizations and surface criteria. 

We now discuss the statistics of the energy levels that 
we find around Ep. In each spectrum, we select a sym- 
metric energy interval (of ~ AOQmeV) around Ep con- 
taining p = 200 levels. Then we find the nearest level 
energy spacing E t — e t +i—e t for each level \t) in the cho- 
sen set, and the mean level spacing 8 — (l/p)St=i ^t- 
To compare our results with the RMT prediction, we 
build a histogram of the calculated spacings using a bin 
size A = 8/ 5. Figure |l|(a) shows such histograms for 
representative realizations of criteria Ci and C 2 . The en- 
ergy is in normalized units: s = E t /8. In RMT our 
system falls into the orthogonal ensemble due to its time 
reversal symmetry. In this case, the Wigner distribution 
P\v(s) — (~Ks/2)exp(— 7rs 2 /4) is a good approximation to 
the distribution predicted for the Gaussian orthogonal 
ensemble (GOE) for the probability of finding a spacing 
E t /8 in the interval (s, s + ds] Q. We see from Fig. ffl(a) 
that the probability of finding a spacing Et/8 in the inter- 
val (s, s + A] (histogram) is well described by P\y(s) xpA 
(thin-solid line) irrespective of the surface criteria. 

In general, given a set of q levels, neighboring level dis- 
tributions Pk(s) are defined as the probability of finding 
a spacing E^ containing k adjacent levels in the interval 
(s,s + ds]. The case k = corresponds to the nearest 
level spacing distribution discussed above. The statisti- 
cal properties of -Pfc(s) are well established for the GOE 
I]. In particular, specific results hold for the standard 
deviation (or width) Ok of these distributions We 
now compare our numerical findings with the GOE pre- 
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dictions. Formally, (a k ) 2 = (<5" 2 /g) J2Li( E t ~ fa) 2 for 
the k-spacing: E\ = Ek+t+i —£t between levels \k + t+l) 
and Sk the mean of Eg, and S = So- Note that ac- 
cording to this definition o~k is expressed in units of S. 
Figure |l|(b) shows the GOE prediction (solid line) for a> 
and the values (symbols) found with the above set of p 
calculated energy levels for different disorder realizations. 
For both surface criteria, the calculated values of at scat- 
ter around the GOE prediction. For criterion C 2 , where 
surface is a large fraction of the particle, the scatter is 
stronger. Although the number of disorder realizations 
shown is too small for accurate ensemble averaging, the 
calculated ak are consistent with the RMT trend. 

We now consider the details of the electronic structure 
in a smaller region around Ef. Figure ^ is a plot of the 
first few levels around \Ep) for five disorder realizations 
for each surface criterion. The energy is measured rel- 
ative to Ef in units of Sg P . The levels present a rich 
structure and are far from being equally-spaced by the 
energy 5g l p . We find that most of the levels form groups 
of two or more. The results are qualitatively similar for Ci 
and C 2 , however, the details around Ep are set by disor- 
der and surface criteria. Note that the predicted number 
of levels (~ 19) for all disorder realizations (with excep- 
tion to #1 and #6) studied here is fairly consistent with 
the number of levels expected from the SP model in the 
considered energy range: [— 10, 10]J^p. It is important to 
note that the clustering seen in the spectra is taking place 
in the presence of disorder in the metal-oxide interface. 
Although we have shown that in a wider range of energy 
around Ep the probability of finding strong bunching of 
levels is small, this does not mean that pairs or larger 
clusters are avoided in a small neighborhood around Ep. 

We now turn to the implications of this clustering of 
energy levels for experimental tunnel resonance spectra 
of the nanoislands. The orthodox model of the device 
(nanoisland+tunncl junctions+leads) treats the leads as 
reservoirs of electrons with chemical potentials [1l and 
fin, coupled capacitively (with capacitances Cl and Cr) 
to the nanoisland through the tunneling junctions. The 
applied voltage V between the leads is divided between 
Cl and Cr setting the the electrochemical potential in 
the right(left) lead to = n R{L) + (-)(C L(fl) /C s ) eV 

with Cs = Cl+Cr. When V = 0, /i£ and /ip are consid- 
ered aligned with Ep. As discussed by Averin and Ko- 
rotkov pi ] , the discrete energy spectrum of the nanois- 
land allows tunneling to occur whenever (measured 
from Ef) matches the energy of one of the available 
levels of the nanoisland. Figure |3| shows the predicted 
flBf dl/dV spectra for different disorder realizations and 
surface criteria. The resonances are marked with flat- 
top(bottom) segments indicating that has origi- 

nated the resonance. Two different regimes: i) (Fig. 
|(a)) symmetric (C L = Cr), and ii) (Fig. |(b)) non- 
symmetric (Cl — l.baF and Cr = 3.2aF ^|) tunneling 
junctions were considered. The dl/dV spectra arising 
from a model with equally spaced energy levels is shown 
(solid circles) for comparison. Clearly, the clustering in 



the energy spectra creeps into the dl/dV spectra in a 
non-trivial way. The nature of the clusters of resonances 
in the dl/dV spectrum is strongly dependent on disorder 
and on the relative values of Cl and Cr. 

Now we show that the clustering of resonaces present 
in the calculated dl/dV qualitatively resembles the ob- 
served clustering. Experimentally (Fig. 1(a) of Ref. ||), 
the clusters of resonances observed are as follows: i) a sin- 
gle resonance at the begining of the spectrum (this sets 
the origin of energies), ii) 4 near-by resonances around 
Sg l p, iii) 2 near-by resonances located at 1.5Sg P , and iv) 4 
resonces distributed over an energy interval equal to Sg l P 
centered around 2.55 g l P . We choose spectrum #4 in the 
symmetric and non-symmetric regime, under criterion Ci, 
for comparison with the experiment. In the symmetric 
regime the spectrum shows: i) a single resonance at the 
origin, ii) 2 pairs of resonances located around Sg l P and 
2(5 g P , respectively, and iii) 4 near-by resonances spanning 
an energy interval approximately equal to Sg l P , around 
3.55g l P . When switching to the non-symmetric regime 
(the experimental regime) the structure of the clustering 
is slightly modified, with the overall number of resonces 
smaller. We see, then, that the observed and calculated 
clusters show qualitative common features. Whether the 
energy scales (cVgp) on which the theoretical and exper- 
imental clusterings occur are similar or not remains an 
open question: This depends on what the correct value of 
the effective volume of the Al island is. We have shown 
that this volume may be significantly larger than has 
been assumed previously. If this is the case then the 
experimental and theoretical energy scales may be very 
similar; letting the experiment be explained in a new way. 

Finally, if we further assume that non-equilibrium 
many-body effects such as those described by Agam et 
al. |j| are present then still more clustering appears in 
the spectrum as most of the single-electron transport res- 
onances in Fig. H split to become clusters of resonances. 
However, our results indicate that the observed tunnel- 
ing resonance spectra should show a different behavior 
from what would follow from the simple equidistant-level 
model of the underlying single-electron spectrum that has 
been assumed in previous analyses of the experiments, 
whether many-body non-equilibrium effects are present 
or not. 

Support from NSERC and CIAR is acknowledged. 
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label 


No 


N Al 
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E F (eV) 


8 (meV) 


#1 


444 


461 


9093 


8.297 


2.34 


#6 


651 


623 


9714 


8.483 


2.32 


#2 


448 


457 


9105 


8.313 


2.30 


#7 


634 


640 


9663 


8.370 


2.30 


#3 


451 


454 


9114 


8.330 


2.30 


#8 


632 


642 


9657 


8.383 


2.32 


#4 


440 


465 


9081 


8.293 


2.31 


#9 


630 


644 


9651 


8.348 


2.34 


#5 


434 


471 


9063 


8.251 


2.32 


#10 


629 


645 


9648 


8.336 


2.35 



TABLE I. Parameters characterizing the nanoparticles. Ci and Ci are criteria that determine the surface (see text). 



FIG. 1. (a) Probability distribution (histograms) for the nearest level spacings and Wigner distribution (thin solid line), (b) 
Standard deviation at for different disorder realizations (points) and the GOE prediction (solid line). Ci and C 2 as in Table 1. 



FIG. 2. Energy levels for different disorder realizations (see text). Ci and C2 as in Table 1. 



FIG. 3. dl/dV spectra for different disorder realizations and capacitance regimes (see text). Ci and C 2 as in Table 1. 
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Figure 1 , Narvaez and Kirczenow 
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Figure 2, Narvaez and Kirczenow 
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Figure 3, Narvaez and Kirczenow 



